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Abstract 

Background: Although sexual reproduction is dominant within eukaryotes, asexual reproduction is widespread and 
has evolved independently as a derived trait in almost all major taxa. How asexuality evolved in sexual organisms is 
unclear. Aphids, such as Acyrthosiphon pisum, alternate between asexual and sexual reproductive means, as the 
production of parthenogenetic viviparous females or sexual oviparous females and males varies in response to 
seasonal photoperiodism. Consequently, sexual and asexual development in aphids can be analyzed simultaneously 
in genetically identical individuals. 

Results: We compared the transcriptomes of aphid embryos in the stages of development during which the 
trajectory of oogenesis is determined for producing sexual or asexual gametes. This study design aimed at 
identifying genes involved in the onset of the divergent mechanisms that result in the sexual or asexual 
phenotype. We detected 33 genes that were differentially transcribed in sexual and asexual embryos. Functional 
annotation by gene ontology (GO) showed a biological signature of oogenesis, cell cycle regulation, epigenetic 
regulation and RNA maturation. In situ hybridizations demonstrated that 16 of the differentially-transcribed genes 
were specifically expressed in germ cells and/or oocytes of asexual and/or sexual ovaries, and therefore may 
contribute to aphid oogenesis. We categorized these 16 genes by their transcription patterns in the two types of 
ovaries; they were: i) expressed during sexual and asexual oogenesis; ii) expressed during sexual and asexual 
oogenesis but with different localizations; or iii) expressed only during sexual or asexual oogenesis. 

Conclusions: Our results show that asexual and sexual oogenesis in aphids share common genetic programs but 
diverge by adapting specificities in their respective gene expression profiles in germ cells and oocytes. 



Background 

Sexual reproduction involves two main events: meiosis 
and fertilization, and creates new genotypes by shuffling 
allelic combinations. Although the predominance of sex- 
ual reproduction in eukaryotes supports this innovation 
as a successful reproduction strategy, asexuality has 
evolved independently multiple times from sexual ances- 
tors in almost all major taxa [1-4], such as in stick 
insects [3] and Ranunculus plants [4]. How asexuality 
has evolved in sexual organisms is unclear. In aphids, 
asexuality was acquired once about 250 million years 
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ago by a common sexual ancestor [5]. Most aphid 
species alternate between sexual reproduction and asex- 
ual parthenogenetic reproduction according to seasonal 
variations. In spring and summer, aphids reproduce 
asexually by parthenogenesis and produce clonal parthe- 
nogenetic female progeny by viviparity. The autumnal 
shortening of the photoperiod induces the concentration 
of juvenile hormone (JH) to decrease in the aphid hae- 
molymph., and particular form of parthenogenetic 
female called the sexuparae are produced. Sexuparae 
females produce sexual females and males that subse- 
quently mate to produce overwintering eggs. Although 
parthenogenetic viviparous females and sexual oviparous 
females exhibit major differences in morphology 
and behavior, they share the same genome. This 
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phenomenon, called reproductive polyphenism, is an 
example of aphid phenotypic plasticity [6]. 

The cellular and cytogenetic bases of reproductive poly- 
phenism have been described for several aphid species 
[7,8]. In aphids, three generations are represented within 
one viviparous female: the mature embryos developing 
inside the maternal abdomen carry the first developmen- 
tal stage of the third generation. This phenomenon is 
known as the 'telescoping of generations' [9]. The 
embryonic developments of asexual and sexual females 
are similar until the formation of undifferentiated ger- 
maria, which takes place at the end of germ band retrac- 
tion (stage 18) [8]. The 32 germ cells within the germaria 
enter into early meiotic prophase (zygotene) and a partial 
association between bivalents occurs [7]. After this stage, 
the embryonic development and maturation of germ cells 
diverge between sexual and asexual females. Ancestral 
sexual development occurs under short photoperiod: the 
basal cells of germaria progress through meiosis; the sex- 
ual oocyte displays chromosome pairing, synapsis and 
meiotic recombinations. Meiosis is blocked in first meio- 
tic prophase until the egg is fertilized and laid. Under 
long photoperiod, sexual oogenesis is inhibited [7]. The 
asexual oocyte ovulates in the posterior part of the ger- 
maria and forms a new follicle where it skips the first 
meiotic division and undergoes a single maturation divi- 
sion resulting in a discarded polar body and a diploid clo- 
nal oocyte. Synchronous mitotic divisions follow 
immediately without fertilization to initiate early embryo- 
genesis, which continues within the parthenogenetic 
female abdomen and progresses through a series of 20 
developmental stages [8]. 

We identified genes regulated in sexual and asexual 
aphid oogenesis by combining transcriptomic analyses 
and in situ hybridizations. In aphids, sexual and asexual 
reproduction can be analyzed in individuals possessing 
identical genomes. Sequencing of the LSR1 pea aphid 
genome [10] and the recent development of new geno- 
mic resources offer an unique opportunity to under- 
stand the molecular bases of sexual and asexual 
reproduction. Our results demonstrate that sexual and 
asexual reproduction pathways require unique and 
shared genetic programs. 

Results 

Synchronization of sexual and asexual embryos 

The development of aphid embryos under long and short 
photoperiods is not synchronous. However, comparing the 
embryonic transcriptomes produced during sexual and 
asexual oogenesis requires synchronized embryo develop- 
ment. Thus, we used kinoprene (a JH analogue) to syn- 
chronize the development of sexual and asexual embryos 
under otherwise identical environmental conditions. Topi- 
cal applications of kinoprene induces sexuparae reared 



under short photoperiods to produce asexual embryos 
instead of sexual embryos [11]. Compared to acetone- 
treatment (control), kinoprene treatment induced com- 
plete reversion of sexual to asexual reproduction in all lar- 
vae laid during the first two days after the onset of 
larviposition (Figure 1). Subsequent transcriptome com- 
parisons were thus restricted to the five most-developed 
embryos dissected from synchronized kinoprene- or acet- 
one-treated sexuparae. These embryos corresponded to 
the progeny laid on the first day. 

Differentially expressed transcripts 

Synchronized asexual (kinoprene-treated) and sexual 
(acetone-treated) embryos corresponding to embryonic 
developmental stages 18, 19 and 20, were collected from 
sexuparae 24 h, 48 h and 72 h after treatment, respectively. 
These are the final three developmental stages in aphid 
embryogenesis and correspond to eye differentiation 
(stage 18), muscle formation (stage 19) and the mature 
embryo (stage 20) [8]. We compared the transcriptomes of 
these embryos using microarrays. Pairwise correlations 
between the five biological replicates were high (r > 0.88), 
indicating good reproducibility. We compiled a list of 33 
transcripts that were statistically differentially expressed: 
none of these transcripts were differentially regulated at 
developmental stage 18. Only one transcript {orb) was 
more abundant in asexual embryos than in sexual embryos 
at developmental stage 19, whereas 33 transcripts (includ- 
ing orb) were differentially-regulated at developmental 
stage 20, including orb. At this final stage, four transcripts 
were more abundant in sexual embryos than in asexual 
embryos, with fold-change (FC) values between 1.49 and 
2.86. Twenty-nine transcripts were more abundant in 
asexual embryos than in sexual embryos (FC values 
between 1.69 and 5.70) (Table 1). 

Based on gene annotation in other organisms, 4 func- 
tional groups were defined (Table 1). Seven transcripts 
were related to oogenesis (e.g. nanos-1, orb, kelch); five 
transcripts were related to post-transcriptional regula- 
tion, such as mRNA polyadenylation (e.g. gld2) or 
mRNA silencing via RNA interference (cbp20); four 
transcripts were related to epigenetic regulation, such as 
histone modifications (suv4-20, uhrfl) or histone var- 
iants (H2B.3, HI); and three transcripts were assigned to 
cell cycle regulation (e.g. cyclin J, clasp 1). Fourteen tran- 
scripts were not assigned to any functional category: 
among them, ten showed no similarity with character- 
ized transcripts. 

Our results displayed multiple genes related to oogesis 
(7/33). This oogenesis signature suggests that our 
experimental design, using artificially-induced asexual 
aphids, successfully targeted genes defining asexual or 
sexual oogenesis within the germaria or the first ovu- 
lated oocytes of the developing embryo. 
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Figure 1 Effect of prenatal application of kinoprene on the type of embryonic ovaries. Acetone or kinoprene dissolved in acetone were 
topically applied to sexuparae individuals, 24 h after the fourth instar molt. Progeny batches were collected after initiation of larviposition, for 
five days every 24 h. The reproductive phenotypes of the progeny were determined and the percentage of each morph (asexual females, 
ambiphasic females that contain a mixture of eggs and embryos, sexual females and males) is indicated. In control conditions (acetone), 
sexuparae sequentially produced sexual females and males: only sexual females were produced during the first three days after the onset of 
larviposition; males appeared the fourth day. Sexuparae treated with kinoprene produced exclusively parthenogenetic progeny the first day after 
the onset of larviposition. On day two, almost all the progeny were parthenogenetic. On day three, the majority of the progeny were still 
parthenogenetic. On days four and five, sexual morphs represented the majority of the progeny. 



Localization of regulated transcripts in sexual and asexual 
ovaries 

We compared the localization of these differentially- 
expressed transcripts in the ovaries of untreated asexual 
and sexual aphids using whole-mount in situ hybridiza- 
tion (WISH), a technique that is not suitable for tran- 
script quantification. Specific antisense riboprobes were 
produced for 25 transcripts: five produced no signal 
(Additional file 1), four showed ubiquitous distribution 
(Additional file 2) and 16 transcripts were specifically 
expressed in germaria (containing germ and nurse cells) 
and/or in oocytes. These results on untreated aphids 
were consistent with the microarray data from acetone- 
or kinoprene-treated aphids: most of the candidate 
genes are expressed in the ovarioles. These 16 tran- 
scripts were grouped based on their patterns of expres- 
sion in the two types of ovaries. 

The first group included nine transcripts (bicC, 
histoneHl, glel, gld2, pop2, arl6ipl, ACYPI007465, 
ACYPI0 10052 and ACYPI25088) that showed the same 
localization in sexual and asexual ovaries (Figure 2) 
despite microarray results that indicated higher levels of 
expression in asexual embryos. For five transcripts (glel, 
bicC, ACYPI2S088, ACYPI007465 and ACYPI010052), 
residual signals were observed in the developing embryo 
germ band (Figure 2). 

The second group included two transcripts (uhrfl and 
orb) that exhibited different localizations in the two 
types of ovaries. Uhrfl signal was detected in germaria 



and oocytes in asexual ovaries, and a weak signal was 
detected in sexual oocytes (Figure 3). orb was detected 
in oocytes prior to ovulation (stage 0), in asexual ger- 
maria, and in unknown cellular structures located at the 
junction between sexual germaria and oocytes. Fluores- 
cent detection of orb riboprobes was used together with 
fluorescent staining of actin filaments and DNA in 
order to define the intracellular localization of orb tran- 
scripts (Figure 4). In asexual ovaries, orb signals could 
be seen during early oogenesis, in stage 0 (oocyte forma- 
tion) and in stage 1 (separation of the oocyte from the 
germarium) (Figure 4A). In sexual ovaries, the tran- 
scripts were detected only before ovulation at the junc- 
tions between the sexual germaria and the oocytes 
(Figure 4A). According to actin staining, a group of cells 
is localized at the junctions between sexual germaria 
and oocytes, these cells may correspond to the pre- 
sumptive oocytes (Figure 4B). Thus, in asexual germaria, 
orb transcripts were detected in oocytes but not in pre- 
sumptive oocytes; whereas in sexual ovaries, these tran- 
scripts were detectable only in cells that probably 
corresponded to presumptive oocytes but were not 
detectable in mature oocytes. 

The third group included five transcripts {Isdl, lode- 
star, cyclin J, ACYPI39770 and ACYPIS4656) that could 
be localized only in one type of ovary (Figure 5). These 
specific expression patterns were confirmed by microar- 
ray expression data for Isdl, lodestar, cyclin] and 
ACYPI39770. The Isdl riboprobe showed a strong 
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Table 1 Description of the 33 differentially-expressed transcripts in sexual and asexual embryos 

Functional group Gene identity Similarity Commentaries 



Fold 
change 



Oogenesis (7) 



Post- 

transcriptionnal 
regulation (5) 



Epigenetic 
regulation (4) 



Cell cycle (3) 



Other (4) 



No similarity (10) 



ACYPI009492-RA nudel 

ACYPI56634- orb 
RA_0 

ACYPI0088 1 1 -RA bicaudal-C 

ACYPI31336-RA kelch 

Ap-nanosl nanos 



ACYPI23247-RA 
ACYPI46386-RA 

ACYPI00753-RA 



putative kelch-like 

putative kelch-like 
repeat 

lodestar 



ACYPI003103-RA Gld2 

ACYPI50630-RA Glei 

ACYPI001842-RA Pop2 

ACYPI49482-RA cbp20 



ACYPI002182- 
RA 

ACYPI005639- 
RA 

ACYPI28709-RA 

ACYPI007975- 
RA 

ACYPI006224-RA 

ACYPI007770-RA 

ACYPI009671-RA 

ACYPI48893-RA 

ACYPI48336-RA 

ACYPI34416-RA 

ACYPI010082-RA 

ACYPI50301-RA 

ACYPI39347-RA 

ACYPI38914-RA 

ACYPI54656-RA 

ACYPI25970-RA 

ACYPI010052-RA 

ACYPI007465-RA 

ACYPI25088-RA 

ACYPI39770-RA 

ACYPI005121- 
RA 



histone H2B.3 

histone H1 

Suv4-20H1 
uhrfl 

cyclin J 
clasp 7 
Drp-1 
Arl6ipl 
Lsd-1 

pol protein 
Six4 

no similarity 
hypothetical protein 
predicted protein 
hypothetical protein 
hypothetical protein 
hypothetical protein 
hypothetical protein 
hypothetical protein 
no similarity 
hypothetical protein 



oocyte dorsal/ventral axis formation, egg activation 

oocyte axis specification, mRNA polyadenylation, germ cell 
development 

mitosis, oogenesis, ovarian follicle cell migration 

female germline ring canal formation, karyosome formation, oogenesis 

female meiosis chromosome segregation, germ cell development, 
oogenesis 

female germline ring canal formation, karyosome formation, oogenesis 
female germline ring canal formation, karyosome formation, oogenesis 

termination of RNA polymerase II transcription, female meiosis 

RNA polyadenylation, histone mRNA catabolic process, mRNA 
processing 

poly(A)+ mRNA export from nucleus 

RNA metabolic process, nuclear-transcribed mRNA poly(A) tail 
shortening 

gene silencing by miRNA, production of siRNA involved in RNA 
interference 

chromatin assembly or disassembly, nucleosome assembly 

chromatin assembly or disassembly, nucleosome assembly 

histone methylation, regulation of gene expression, epigenetic 
DNA repair, cell cycle, methylated histone residues binding 

mitotic cell cycle, embryonic 

cell division, exit from mitosis, condensed chromosome kinetochore 

cytokinesis, intracellular distribution of mitochondria 

cotranslational protein targeting to membrane 

lipid particle organization, regulation of lipid storage 

RNA dependant DNA Polymerase/Reverse Transcriptase 

fat body development, gonad development, mesoderm development 



3.86 
4.67 

4.69 
1.81 
3.69 

1.69 
-1.49 

3.39 

4.03 

2.46 
2.21 

2.11 

3.60 

2.85 

1.75 
2.48 

5.70 
1.48 
1.79 
4.25 
-2.86 
-1.92 
1.56 
3.92 
2.70 
2.04 
2.12 
3.86 
3.78 
2.53 
2.20 
2.88 
-1.78 



Positive fold-change values indicate transcripts that were more abundant in asexual embryos. Negative fold-change values indicate transcripts that were more 
abundant in sexual embryos. 



specific signal in the germaria and oocytes of sexual 
ovaries. In contrast, a specific signal for lodestar, cyclin 
J, ACYPI39770 and ACYPI54656 was detected only in 
asexual ovaries (Figure 5). lodestar transcripts were 
restricted to the diploid oocyte after separation from the 



germarium (stage 1). cyclin J, ACYPI39770 and 
ACYPIS46S6 transcripts were detected in asexual ger- 
maria and oocytes (Figure 5). In addition, ACYPI54656 
and ACYPI39770 exhibited a residual signal in the germ 
bands of developing embryos (Figure 5). 
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Figure 2 Transcripts with identical localizations in both sexual and asexual ovaries. Histone HI, glel, gld2, bicC, pop2, arl6ip1, ACYP125088, 
ACYPI007465 and ACYPI10052 transcripts were detected in germaria (hollow arrowheads) and oocytes (black arrowheads) of both sexual and 
parthenogenetic ovaries. Sense riboprobes were used as negative controls. For glel, bicC, ACYPI25088, ACYPI007465 and ACYPI010052, a residual 
signal was observed in the developing embryo germ band (arrows). Bar scale: 50 urn. 
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Figure 3 Transcripts with different localizations in sexual and asexual ovaries. Uhrfl transcripts were detected in germaria (hollow 
arrowheads) and oocytes (black arrowheads) of asexual ovaries whereas a weak signal was detected in sexual oocytes (black arrowheads). Orb 
transcripts were detected in the basal parts of germaria in asexual and sexual ovaries (hollow arrowheads) with distinct localizations. In asexual 
germaria, these cells correspond to oocytes, whereas the cellular structures containing orb transcripts in sexual ovaries are undetermined. Sense 
riboprobes were used as negative controls. Bar scale: 50 urn. 
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Figure 4 Detection of orb transcripts in sexual and asexual ovaries by fluorescent WISH. (A) orb riboprobes were detected with an 
antibody coupled to alexin fluorochrome (red) and DNA was stained with DAPI (white). In sexual ovaries (upper pictures), orb transcripts were 
detected at the basal part of the germaria (hollow arrowheads). In asexual ovaries (lower pictures), orb transcripts were detected in oocytes 
(white arrowheads) before ovulation (stage 0), and after ovulation (stage 1). (B) Actin filaments and DNA in sexual ovaries were stained with 
phalloidin (green) and propidium iodide (red), respectively. Cells at the basal part of the sexual germaria possibly correspond to presumptive 
sexual oocytes (double arrows). Orb sense riboprobes were used as negative controls. Bar scale: 30 urn. 



Discussion 

By coupling quantitative transcriptomic and qualitative 
WISH approaches, we identified 33 transcripts asso- 
ciated with the onset of divergent transcripts expression 
leading to the establishment of the sexual or asexual 
phenotype in aphids. The involvement of these candi- 
date genes in determining the developmental direction 
of oogenesis is supported by: i) a strong signature of 
oogenesis among the regulated transcripts (7/33); and ii) 
the oocyte- and/or germarial-specific localizations of 16 
of the transcripts. These data illustrate both shared and 
unique regulatory patterns of transcript expression for 
the two modes of reproduction. 

First, we identified transcript expression profiles 
shared by sexual and asexual oogenesis. These 



transcripts are expressed in sexual and asexual germaria 
and/or oocytes but display differences in their expres- 
sion levels and/or in their localization. In various organ- 
isms, gld2, bicC, pop2 and orb are all involved in 
regulating the poly(A) tail length of maternal mRNAs. 
This process is a common mechanism of translational 
control essential for progression of meiosis and axis pat- 
terning during oocyte maturation [12-14]. The regula- 
tion of poly(A) tail length of maternal mRNAs results 
from a balance between concomitant deadenylation 
(translation repression) and polyadenylation (translation 
activation) [15]. In Drosophila melanogaster, gld2, bicC 
and orb have been described as part of the cytoplasmic 
polyadenylation complex that polyadenylates mRNAs 
and activates their translation [13,16]. Gld2-dependent 
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Figure 5 Transcripts detected specifically in one type of ovary. Lsdl transcripts were detected specifically in the germaria (hollow 
arrowheads) and haploid oocytes (black arrowheads) of sexual ovaries, lodestar, cyclin J, ACYPI39770, ACYPI54656 transcripts were detected 
specifically in asexual ovaries. Lodestar signal was detected specifically in asexual oocytes (black arrowheads). Cyclin J, ACYPI54656 and ACYPI39770 
signals were detected in germaria (hollow arrowheads) and in oocytes (black arrowheads) of asexual ovaries. ACYPI54656 and ACYPI39770 showed 
a residual signal in the germ bands (arrows) of developing embryos. Sense riboprobes were used as negative controls. Bar scale: 50 urn. 



mRNA polyadenylation and translational activation have 
been shown to be essential for oogenic meiosis [13]. On 
the other hand, mRNA deadenylation by the major 
deadenylation complex in Drosophila, CCR4-NOT, 
which includes the deadenylases CCR4, POP2, and four 
NOT proteins [17], promotes negative regulation of tar- 
get mRNAs. Deadenylation of specific maternal mRNAs 
could be involved in the precise temporal and spatial 
mRNA localization necessary for axial patterning during 



oogenesis [12]. Expression of Api-gld2, Api-bicC, Api- 
pop2 and Api-orb in sexual and asexual germaria and/or 
in oocytes suggests that the balance between cytoplas- 
mic polyadenylation and deadenylation in oocytes acts 
in the maturation of haploid as well as diploid oocytes 
in the pea aphid as it has been shown for other sexual 
organisms. Moreover, expression of the oocyte-specific 
linkers histoneHl and uhrfl in asexual and sexual 
oocytes suggests that epigenetic regulation is involved in 
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pea aphid oogenesis. The linker histone HI interacts 
with DNA to set-up and maintain the organization of 
the chromatin between the nucleosomes. Histone HI 
proteins form a complex family of related proteins that 
have distinct species, tissue and developmental specifici- 
ties and are involved in the epigenetic control of gene 
expression. The replacement of a somatic linker HI by 
an oocyte specific one (Hloo) during oogenesis has been 
described in a wide range of species (reviewed in 
[18,19]) and may promote a loosening in chromatin 
structure that appears to correlate with the initiation of 
meiosis [19]. The presence of Api-histoneHl in germaria 
and oocytes suggests that this oocyte-specific variant of 
HI might be involved in both sexual and asexual oogen- 
esis in aphid, uhrfl encodes an ubiquitin-like protein 
that is involved in two epigenetic silencing pathways 
through its regulation of DNA and histone methylation 
[20,21]. Although, uhrfl has not been implicated in any 
oogenesis process to date, our results suggest a role for 
this protein in asexual and sexual oogenesis in aphids. 
Altogether, our results suggest that epigenetic regulatory 
mechanisms and the balance between cytoplasmic 
mRNA polyadenylation and deadenylation may have 
roles in regulating haploid and diploid oocyte matura- 
tion. Moreover, as orb transcript is the earliest to be dif- 
ferentially expressed during the developmental time 
course, we propose that Orb might play a major role in 
determining oocyte fate. 

Second, we identified one gene, Isdl, that is specifi- 
cally expressed in pea aphid sexual germaria and 
oocytes. The Isdl transcript was not significantly 
detected either in asexual embryos according to micro- 
array data (Additional file 3) or in asexual ovaries 
according to WISH. Isdl encodes a protein involved in 
the regulation of lipid droplet storage in eukaryotic cells 
[22]. Lipid droplets, like mature yolk bodies, are the 
major energy storage for oocytes and developing 
embryos in Drosophila [23]. Our results suggest that 
Isdl is required only for sexual oogenesis during which 
it may be involved in yolk accumulation. The absence of 
Isdl transcripts in asexual embryos is consistent with 
the absence of yolk in parthenogenetic oocytes [8]. 
Thus, identification of Isdl expression as a specific fea- 
ture of sexual oogenesis validates the accuracy of our 
approach. Despite the existence of genes known to be 
involved in meiotic recombination in the pea aphid gen- 
ome [24], we identified no such genes as specific mar- 
kers of sexual oocytes. However, germline expression of 
meiotic genes such as spoil was shown to be similar in 
sexual and asexual aphids. Expression of this gene 
seems to be modulated between asexual and sexual 
aphids by alternative splicing [25]. 

Finally, we identified three transcripts that were specific 
to asexual oogenesis. No significant transcription signal 



has been observed for these genes either in sexual 
embryos using microarrays (Additional file 3) or in sexual 
ovaries using WISH. That the polypeptide encoded by 
ACYPI39770 does not bear any similarity to a protein or 
protein domain of any other organism, suggests that 
ACYPI39770 is either an orphan gene or a gene with an 
unusual accelerated evolution rate. In Drosophila, cyclin J 
and lodestar are involved in sexual oogenesis and early 
embryogenesis [26,27]. In aphids, their expression profiles 
suggest that they may haveroles in asexual oogenesis 
only. The absence of lodestar and cyclin J transcripts 
from aphid sexual oocytes and embryos was particularly 
unexpected. Three lodestar paralogs were identified within 
the pea aphid genome {ACYPI000753, ACYPI000590 
and ACYP1006971). Although the three copies were 
represented on the microarray by copy-specific oligonu- 
cleotides, only ACYPI0007S3 was differentially expre- 
ssed in sexual and asexual embryos. Localization of 
ACYPI0007S3 transcripts by in situ hybridization showed 
that its expression was specific to asexual oocytes. How- 
ever, although the sequence of the riboprobe used to 
locate ACYPI0007S3 transcripts was absent from the 
ACYP1006971 sequence, it shared 87% nucleotide identity 
with the unregulated paralog ACYP1000590. Therefore, 
cross hybridization between the ACYP1000753 riboprobe 
and ACYPI000S90 transcripts cannot be excluded. Never- 
theless, we detected no ACYPI0007S3 expression in sex- 
ual oogenesis. This unexpected result might reflect the 
functional specialization of at least one of the lodestar 
paralogs in asexual oogenesis. 

Among the 16 transcripts found to be specifically 
localized in germaria and/or oocytes , Isdl was the only 
one for which the expression level was found to be 
higher in sexual ovaries than in asexual ovaries accord- 
ing to the microarray results. This disproportion is con- 
sistent with the striking over-representation of 
transcripts that were upregulated in asexual embryos 
(29/33), as evidenced by the microarray analysis. Sexual 
oogenesis and asexual oogenesis display major inherent 
differences. Within embryos, the sexual oocyte nucleus 
remains blocked in meiosis prophase 1 whereas the 
asexual oocyte undergoes a modified but complete 
meiotic division immediately after ovulation [7]. This 
modified meiosis involves specific changes within the 
oocyte; self-organized asters are formed, that then 
recruit centriole precursors and pericentriolar material 
[28]. This process may be accompanied by the activation 
of microtubule stabilization factors or by the inactiva- 
tion of microtubule destabilization factors specific to 
asexual oogenesis [28]. The predominance of genes that 
are upregulated in asexual embryos may be partially 
explained by the production of molecular components 
that are specifically required for asexual oogenesis but 
not for sexual oogenesis. 
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Conclusions 

In this study, we identified a number of putative embryo- 
nic target genes of the signaling cascade initiated by 
photoperiodic cues and transduced by JH. These genes 
are likely to be among those that direct sexual or asexual 
differentiation of the developing embryos. Our analysis of 
their expression patterns represents the first time the 
molecular bases of aphid sexual and asexual oogenesis 
have been addressed. Altogether, our results showed that 
aphid asexual oogenesis requires the same genes as sex- 
ual oogenesis but that these two distinctive differentia- 
tion processes display specific gene expression programs. 

Methods 

Aphid rearing, dissection and extraction 

Clonal descendants from the sequenced Acyrthosiphon 
pisum LSR1 clone were reared on Vicia fabae at 18°C 
and at low density (5 individuals per plant) to prevent 
the production of winged morphs. Parthenogenesis was 
maintained under long photoperiod (16 h). In our con- 
ditions, the shortening of the photoperiod (12 h) 
requires three generations to induce the production of 
sexual morphs. The first asexual generation was reared 
under long photoperiod and transferred to short photo- 
period (12 h) after the third instar moult. Adult asexual 
virginoparae females produced the second asexual gen- 
eration, named sexuparae. Adult sexuparae, reared in 
short photoperiods, produced the third generation of 
sexual females oviparous, males, and asexual viviparous 
females. 

Kinoprene application 

Aphid sexuparae were synchronised at the fourth instar 
moult, during a 6 hour window. 24 h after, 400 ng of 
kinoprene (Sigma Aldrich) diluted in 50 [iL of acetone 
were ectopically applied on sexuparae abdomens of 100 
females. The progeny of treated aphids was collected 
daily for 5 days. When they reached adulthood, collected 
aphids were dissected and their reproductive type was 
determined according to the type of ovaries they con- 
tained within their abdomen: sexual ovaries that are 
only made of haploid eggs, ambiphasic ovaries that con- 
tain a mixture of eggs and embryos, and asexual ovaries 
that are only constituted of embryos at different stages 
of development. Ectopic applications of 50 \iL of acetone 
were performed as a negative control. 

Sampling, RNA isolation and microarray hybridization 

Synchronized sexuparae were randomly divided into 2 
batches and treated with kinoprene or acetone (as a 
control) 24 h after fourth instar moult (Additional file 
4). 25 treated sexuparae were collected 24, 48, and 72 
hours after kinoprene or acetone application. For 
each condition, the 5 most developed embryos were 



isolated from each of the 25 treated sexuparae by dis- 
section, pooled together, frozen into liquid nitrogen 
and stored at -80°C. These most developed embryos 
within treated sexuparae collected 24, 48, and 72 
hours after treatment correspond to the develope- 
mental stage 18, 19 or 20 respectively. This procedure 
was repeated 5 times to generate as much indepen- 
dent biological replicates. Total RNAs were isolated 
from each sample by using the RNeasy Mini kit (Qia- 
gen) according to manufacturer's instructions. RNA 
quality was checked on Bioanalyser (Agilent) and 
quantified on Nanodrop (Thermo scientific). For each 
sample, 20 [ig of total RNAs were sent to the Nimble- 
Gen expression array platform (Roche). Double 
stranded cDNA synthesis and Cy3 end-labelling were 
performed by NimbleGen. 

Microarrays design and analysis 

Custom microarrays were constructed on NimbleGen 
(Roche) 385K 4-plex (4 x 72,000 probes). 24,011 tran- 
scripts were represented by 3 60-mers oligonucleotides 
probes (MIAMExpress, http://www.ebi.ac.uk/miamex- 
press/, microarray: INRA-BF2I_A.pisum_Nimblegen- 
ACYPI_4x72k_vl; Array express accession: A-MEXP- 
1999). Hybridization and scanning were performed by 
NimbleGen, providing the final raw data file (experi- 
ment: Aphid_embryo, ArrayExpress accession: E-MEXP- 
3481). The limma 2.16 package [29] in R 2.9.2 [30] was 
used for statistical analyses. For each transcript the aver- 
age of the 3 probe signals was considered. Median nor- 
malization between microarrays was performed. A 
design matrix incorporating effects of treatment K 
(kinoprene) or A (acetone) and developing stage (18, 19, 
20) was constructed. Variance was adjusted for Bayesian 
fitting of the model. Differential expression was deter- 
mined with a two-way ANOVA considering the effect of 
time and treatment between 3 contrasts: 'A18 versus 
K18', 'A19 versus K19' and 'A20 versus K20\ Significance 
of differential expression was assigned with a 10% false 
discovery rate (FDR) [31]. Annotation of differentially 
expressed transcripts was curated using pea aphid RNA- 
seq data (published in SRA of NCBI), the software 
CAP3 [32] and Blastx on NCBI non-redundant protein 
sequence data bank. Functional annotations were per- 
formed using the Gene Ontology (http://www.geneon- 
tology.org/). 

Riboprobe synthesis for in situ hybridization 

Templates for synthesis of riboprobes were obtained 
from full-length cDNAs collection (ACYPI003103, 
ACYPI010052 [33]) or amplified by RT-PCR and cloned 
(Additional file 5). Total RNAs were extracted from 
parthenogenetic virginoparae females with RNeasy plant 
kit (Qiagen). DNA contaminations were removed by a 
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treatment with RQ1 RNase-free DNAse (Promega). First 
strand cDNAs were produced from 1 [ig of total RNAs 
by using random primer 9 (New England BioLabs) and 
Superscript® III Reverse Transcriptase (Invitrogen) fol- 
lowing the supplier's instructions. cDNAs were used as 
a matrix for PCR amplification with specific primers 
(Additional file 5). Amplified fragments were cloned 
into the StrataClone PCR Cloning Vector pSC-A-amp/ 
kan (StrataClone) or pENTR Directional TOPO (Invitro- 
gen) and sequenced (Genoscreen). Linear PCR products 
were amplified from cloned sequence with universal pri- 
mers M13 and used as a matrix for synthesis of sense 
and antisense riboprobes by using digoxigenin-labelled 
dNTPs and the appropriate RNA polymerase (T3/T7/ 
SP6) supplied in the DIG RNA labelling kit (Roche). 
Remaining DNA was removed with RQ1 RNase-free 
DNAse treatment (Promega) and riboprobes were puri- 
fied with the RNeasy mini kit (Qiagen). Riboprobe qual- 
ity and quantity were checked on Nanodrop (Thermo 
scientific). 

Whole mount in situ hybridization and microscopy 

Whole mount in situ hybridization was performed on 
ovaries of sexual or asexual nymphs. Ovaries were dis- 
sected and fixed in 4% paraformaldehyde in lx PBS at 
room temperature (RT) for 30 min. Ovaries were 
washed in 50% methanol for 30 min, then dehydrated 
and stored in methanol at -20°C. They were rehydrated 
in graded methanol/PTw (lx PBS, 0.2% Tween-20) solu- 
tions (70%; 50%; 30%, 10 min each), post fixated in 4% 
paraformaldehyde in lx PBS for 20 min and washed 3 
times for 5 min in PTw. Ovaries were washed for 45 
min in 1% SDS, 0.5% Tween-20, 50 mM Tris-HCl (pH 
7,5), 1 mM EDTA (pH 8), 150 mM NaCl and 5 times 
for 5 min with PTw. Pre-hybridization was performed 
10 min in hyb-wash solution (5x SSC, 50% formamide, 
0.1% Tween-20) and 1 hour at 65°C in hyb solution 
(0.3% SDS, 5x SSC, 50% formamide, 100 ug/mL heparin, 
0.1% Tween-20, 100 [ig/mL yeast RNA, 10 mM DTT). 
Hybridization was performed overnight at 65°C with 500 
ng/mL of denatured sense or antisense RNA probe 
diluted in hyb solution. Nonspecific hybridizations were 
washed off twice at 65°C in hyb-wash solution for 30 
min, 3 times in 50% hyb-wash solution/ PTw 30 min and 
3 times in PTw for 10 min at RT. Ovaries were incu- 
bated in blocking solution (0.2% BSA in PTw) for 5 min 
and 1 h at RT, and bound digoxigenin labelled probes 
were detected overnight at 4°C with anti-DIG-alkaline 
phosphatase (AP) Fab fragments (Roche) diluted 1:2000 
in blocking solution. Ovaries were washed 4 times 20 
min and twice for 1 h in blocking solution at RT. After 
3 washes in AP reaction buffer (100 mM Tris (pH 9,5), 
100 mM NaCl, 5 mM MgCl 2 , 0.2% Tween-20), signal 
was revealed with 4 \A NitroBlue Tetrazolium/5-Bromo- 



4-Chloro-3-Indolyl Phosphate (NBT/BCIP) Stock Solu- 
tion (Roche) /ml AP reaction buffer. Finally, ovaries were 
rinsed at least 3 times for 5 min in PTw, dehydrated 5 
min in methanol, and rinsed twice again 5 min in PTw, 
before mounting in 70% glycerol in PBS. Samples were 
photographed with a microscope Nikon 90i connected 
to a Nikon type DS-Ril camera or with Olympus BX61 
connected to Nikon DS-Fil camera. 

For the fluorescence protocol, the same procedure 
was followed until hybridization step. Then TSA™ Bio- 
tin System (Perkin Elmer, Waltham, USA) was used for 
detection as follows. After overnight hybridization of 
the probes, nonspecific probes were washed off twice 
at 65°C in hyb-wash solution for 30 min, once in 50% 
hyb-wash solution/TNT buffer (0.1 M Tris-HCl (pH 
7.5), 0.15 M NaCl, 0.05% Tween-20) for 30 min at RT, 
and 3 times in TNT buffer for 10 min at RT. Ovaries 
were blocked with the TSA™ blocking reagent, for 30 
min at RT. Bound digoxigenin labelled probes were 
detected with HRP conjugated anti-DIG (Boehringer- 
Mannheim) (1:250) diluted in TNT buffer at 4°C for 
overnight incubation. Ovaries were washed four times 
in TNT buffer for 15 min. The TSA™ Biotin System 
Amplification was used by incubation 20 min into a 
diluted solution of the Biotinyl Tyramide (Amplifica- 
tion Reagent) in lx Amplification Diluent (1:50). HRP 
catalyzes the formation of TSA free radicals, which 
form covalent bounds to tyrosine residues proximal to 
HRP. They were then washed 3 times 15 min in TNT 
buffer, and 4 times 20 min in PTw. Ovaries were incu- 
bated 2 h into Alexa 594-conjugated Streptavidin (Invi- 
trogen) diluted in TNT buffer (1:500) and rinsed into 
PTw before performing a nuclear stain with TO- 
PR03® (1:1000 in PTw) for 30 min. Finally, ovaries 
were washed 3 times in PTw before mounting. Images 
were acquired and processed using an Olympus 
FV1000 confocal microscope. 

Phalloidin and propidium iodide staining 

Oviparous ovaries were dissected and fixed as previously 
described. Samples were washed 3 times in PTw, fol- 
lowed by PTw for 30 min and 2% Normal Goat Serum 
(NGS)/BPTw for 60 min. Samples were then incubated 
overnight at 4°C with Alexa Fluor 488 phalloidin (Invi- 
trogen) at a concentration of 1:200 to visualize F-actin. 
The next day samples were washed 4 times in PTw, and 
placed in 20 [ig/mL RNaseA (QIAGEN)/PTw for 3 h at 
RT. Samples were washed twice in PTw, and incubated 
with 10 |ig/mL propidium iodide (Invitrogen) in 20 \igl 
mL RNaseA/PTw for 60 min at RT to visualize DNA. 
Samples were washed 4 times in PTw and mounted in 
Vectashield (Vector Laboratories). Images were acquired 
and processed using an Olympus FV1000 confocal 
microscope. 
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